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ABSTRACT 

In this work Gamma Ray Burst (GRB) data is used to place constraints on a putative 
coupling between dark energy and dark matter. Type la supernovae (SNe la) con- 
straints from the Sloan Digital Sky Survey II (SDSS-II) first-year results, the cosmic 
microwave background radiation (CMBR) shift parameter from WMAP seven year 
results and the baryon acoustic oscillation (BAO) peak from the Sloan Digital Sky 
Survey (SDSS) are also discussed. The prospects for the field are assessed, as more 
GRB events become available. 
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1 INTRODUCTION 

The nature of dark energy and dark matter remains an 
outstanding open problem in cosmology. In spite of the 
success of the ACDM parameterization, one must con- 
sider more complex models in order to cast some light on 
the substance of the dark components of the universe. In 
this work one considers models with interacting dark en- 
ergy and dark matter components. There are several useful 
tools to prob e the phenomenology of thes e models, such as 
CMBR data faomatsu et al.ll201ol ). BA O jReid et al.ll200 ' 
iMoIdenhauer fc Ishakll200g| ). SNe data (jKessIer et alJl200: 
and the deviation from the virial equilibrium of galaxy clus 
ters ilBertolami et a l. 2007. 2009; Abdalla 2009h. It has been 
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suggested (jSchaeferl 120031 : iBertolami fc Tavared 120061 ) that 
GRB may be used to extend the Hubble diagram to high 
redshifts, greater than z = 5. At these epochs the Universe 
was dominated by dark matter, from which follows that this 
tool is less sensitive to dark energy. How ever, for models 
wher e dark energy and matt er are coupled (lAmendola. et al.. 



2003; 



2001; 



Bertolami et al.l 20071) or unified ( Kamcns hcik et al.l 



Bento et al.ll2002l . 120031 ; iBarreiro et alj|200gft. GRBs 
be a particularly usefull tool ( Bertolami fc Tavaresl 



migh t 
l2006h . 

In the late 1960s, the Vela array of military satellites 
detected flashes of radiation originating in apparently ran- 
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dom directions in space. The observed bursts lasted between 
tens of milisecond and thousands of seconds, and were com- 
posed of soft (0.01 to 1 MeV) gamma rays. Subsequently 
space missions such as the US Apollo program and the So- 
viet Venera probes confirmed the existence of the GRBs, 
even though its rate of occurrence was virtually unknown 
until the deployment of the Compton Gamma Ray Obser- 
vatory, in 1991. This observatory was equipped with a sen- 
sitive gamma-ray detector, the Burst and Transient Source 
Explorer (BATSE) instrument which was able to detect one 
or two events per day. The collected data allowed to di- 
vide GRBs into two categories: short duration bursts (short 
bursts) and long duration bursts (long bursts). The former 
usually last for less than two seconds and are dominated by 
high energy photons; the latter last longer than two seconds 
and are dominated by lower energy photons. However, this 
distinction is not always clear. 

The physical origin of GRBs has been debated for a 
long time, before their exact position an d a reliable esti - 
mate of their distance was lacking (see e.g. (|Bcrtolamil ll999h 
and references therein). In 1997, several GRBs were detected 
by the BeppoSAX sattelite. A GRB prompt emission is fol- 
lowed by an afterglow emission composed by all wavelengths. 
Depending on its brightness, an afterglow can last from 
days to months after the burst itself, the transient phase. 
The detection of the afterglow did manifold the information 
on GRBs. Through their afterglow, GRB's X-ray, optical 
and ra dio counterparts w ere observed, as well as their red- 
shifts (|Fan fc Piran|[2008l ). confirming the cosmological ori- 
gin of most, if not all, the GRBs. When, in 2003, the long 
GRB 03032 9 was discovered and linked with the supernova 
SN2003dh (|Cobb et al.l 120041 ). it became clear that GRBs 
are linked with the release of gravitational energy during the 
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collapse of stellar mass objects. GRBs are most likely colli- 



mated, considering they reach integrated luminosities w-fee-j— 
L ~ 10 53 erg s -1 , making it hard to associate them w 
an astrophysical object otherwise. This high energy release 
creates an outflow that expands relativistically. Two forms 
of shocks are ensued by the burst, the forward shock and the 
reverse one, wh ich one separated by a contact discontinuity 
(|Lvutikovll2009l ). If the ejected plasma is too strongly mag- 
netized, only the forward shock is formed. One suggested 
possibility is that the prompt emission is generated in a 
baryon dominated ejecta through internal shocks, while the 
forward and reverse shoc ks yield the lon g lasting broadband 
emission, the afterglow (|Lvutikovll2009l ). Actually, the full 
understanding of the prompt emission mechanism, a basic 
GRB property, is still lacking. One possibility is that the 
prom pt emission consists of synchrotron radiation (|Lvutikovl 
120091 ). by the relativistic charged particles moving on the 
magnetized ejected plasma. Currently, the GL AST/FERMI 
mission, in operation, is continuously increasing the avail- 
able GRB data and making it worth, as will be discussed 
and pursued in this work, considering future prospects for 
the subject. For an ov erview of most recent missions see e.g. 
(|McBreen et alj|2010h and references therein. 

GRBs can be use d as distance in di cators 
jGhirlanda et all 120041: Bertolam i fc Tavaresl l200d 
iLiang fc Zhand l200rj ; lAmati et al.l |2008h . Rs main at- 
tractiveness is that the redshift range extends much higher 
than that of SNe la. The main observables that can be 
measured when studying GRBs are its spherical equivalent 
energy, its peak isotropic luminosity, the peak energy 
of its spectrum, the photon fluence, the energy fluence, 
the pulse duration and the redshift of its host galaxy. 
Several empirical correlations among these variables can be 
established. However, there are still large uncertainties in 
their calibration. Furthermore, there is still no satisfying 
physical mechanism accounting for them, so that assuming 
that they hold true can introduce systematic uncertainties 
in our distance indicator. From the existing correlations, 
the very discussed Ghirlanda relation uses the peak energy 
of t he spectrum,, E v jand t he collimation corrected energy, 
E-f (|Ghirlanda et all 120041 ). On the other hand, the Liang 
and Zang relation correlates the isotropic equivalent energy, 
Ei so wi th _E p .i and the j et br eak time of the afterglow of the 
burst (|Liang fc Zhand l200fj| ) . Finally, the A mati relation, 
corre lates the isotropic energy, Ei so , with _E Pi i (|Amati et al.l 
120081 ). This relation is particularly interesting since the 
E Pt i - Ei so correlation requires only two parameters that can 
be directly inferred from the obervations. This correlation 
further emphasizes the relevance of the GRB data. Notice 
that the aforementioned synchrotron process reproduces 
the Amati correlation, a quite interesting feature. 

In this work, GRB data and the Amati relation, in par- 
ticular, are used to probe a generic dark energy - dark mat- 
ter interacting model. In section [21 the interacting model is 
presented. The Amati E p j - Ei an correlation is introduced 
and discussed in section [3] The set of real GRB data is then 



extended to a mock sample of 500 GRBs using a method 
detailed in subsection 13.11 In section 31 one discusses the 
constrains obtained from SNe data in subsection 14.11 BAO 
in subsection 14.21 and CMBR shift parameter in subsection 
14.31 In section [5] one presents the obtained results. In section 
[6] conclusions are presented. 



2 DARK ENERGY AND DARK MATTER 
INTERACTION 

The cosmological model consists of homogeneous matter 
(dark matter and baryons) and dark energy, where the dark 
matter and energy are interacting and have equations of 
state, pdm = and pde = wpDE, respectively. The cou- 
pled energy densities with a coupling £ evolve as follows 
(|Bertolami et alj|2007h : 



f>DM + 3HpDM = (HpDM , 
pDE + 3HpDE{l + w) = -(HpDM 



(1) 

(2) 



The analysis a ssumes for the ratio of the dark compo- 
nents that (|Bertolami et al.ll2007h 
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(3) 



for a constant 77, where, a, is the scale factor, assumed that 
at present an = 1, and 2 is the redshift. 

ratio Jpiancni 

Inserting the time derivative of Eq. (|3]) into Eqs. ([1]) and 
©, one obtains for the coupling 



Co 



DE„ 



(4) 



where Co = — (ti + 3w)Q,de ■ Note that when r\ = — 3w, there 
is no interaction between dark energy and dark matter, since 

C = o. _ Mai 

The solutions for Eqs. ((1} and ([2]) can be written as 



PDM = (!+«) PDM 



D.DE {l + z) V +n D M 



DE Q 



Q 



DMq 



(5) 



PDE 



(1 



_„ + 3 \Q,de q (X + z ) v + Qdm 

pDE 



Inserting Eqs. ([3]), $5$ and $6$ 



£Ide + Qdm 

IDE 



EM J aplDE ( 6 ) 
and ((6|) into the Friedmann equa- 
tion for a fiat universe, H 2 — Hq(pdm + PDE + pb)/ po, where 
Ho and po are the Hubble constant and the total energy den- 
sity at present, one obtains 



E\z) = (l + z) 3 



M 



(O,DM + &DE ) 



(7) 

where fij, and Qdm are the baryon and DM energy densi- 
ties at present and E(z) = H{z)/ Ho- 
lt is interesting to point o ut that the Gene r alized 
Chaplygin Gas model (GCG) (jKamenshcik et all l200ll : 
iBento et alll2002h . an unified model of dark energy and dark 
matter, can be seen as a particular cas e of this interactin g 
model for 77 = 3(1 + a) and w = -1 (jBento et all 120041 ). 
q being the GCG equation of state parameter, p = — s, 
where A is a positive constant. 



3 GAMMA RAY BURTS 

In this work one considers the Amati E Pl i - Ei so correla- 
tion (for a discussion of the use of other correlations see f or 
instance ijBertolami fc Tavarej2006l : lLiang fc Zhandl2006t )). 
This correlation can be used to place constrains on the Hub- 
ble diagram. The sample provided in (|Amati et al.l 120081 . 



DMDEratio 



etaO 



solDM 



solDE 
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that includes the observations of 95 GRB with mea- 
surements for -Ep.i, -Eiso and redshift, is adopted. 

The value of E P) j is an observable quantity, indepen- 
dent of a cosmological model. On the other hand, E- lso is 
computed for each GRB from its spectral parameters, flu- 
ence and redshift using a specific cosmolog ical model (The 
-Eiso presented in (jAmati et a l. 2008, 20Q3) is computed in 
the context of the ACDM scenario with h = 0.7, SIm = 
Qdm + fit = 0.3 and Hde = 0.7). The Amati corre- 
lation assumes a power law relationship between E Pi i and 
iJiso- A cosmological model can then be tested comparing 
the _Ei so computed from the observations with a theoretical 
Ei so obtained from -E p .i. In practice, however, since the E p ,i - 
Eiso relationship is not calibrated, one has to simultaneously 
fit for the power law paramet ers. 

Furthermore, following (Amati 2006), the scatter of 
the E p ,i-Eiso relation that cannot be explained by statis- 
tical fluctuations alone must b e taken into account. As in 
(|D'Agostinill2005l ; lAmatill2006T) . this is done by introducing 
a third parameter, an extrinsic variance, cr ox t, in the fitting 
of the data. Using a power law E p ,i-Ei so relation, 



log _E p ,i = m log Eiso + q , 
one aims to minimize the likelihood function, where 

(log E p ,i - m log E iso - q) 2 



(8) 

- 



2 

Xgrb 



E 

GRB 

objects 



log (27T<7 



(9) 



ffl*t + °"p + Tn 2 af sn , where 



and the total variance is a 2 
<7 P and Uiso are the observational variances on log _E Pi i an 
log .Eiso, respectively. The fit is performed for the three pa- 
rameters m, q and o" e xt- 

This minimization is then carried out for different values 
of the cosmological parameters, resulting in a profile of the 
likelihood function. 



3.1 Generating a GRB mock sample 

Given that the GRB data is currently rather limited, the 
analysis is extended to include a mock sample of GRBs in 
order to test the efficiency of the Amati relation on con- 
straining dark energy and dark matter interacting models. 
The goal is to check the effect of a larger number of GRBs, 
but also the effect of higher redshift GRBs. From 2009 on- 
wards, most of the useful events to fit the Amati re- 
lation came from the Swift and Fermi experiments. 
One can expect from these experiments approxi- 
mately 10 useful GRB events per year. The future 
launch of the EXIST mission, scheduled for 2017, 
will considerably improve this rate, hopefully simul- 
taneously reducing the measurement errors. This 
paper settles on a best case scenario of 500 GRBs 
events but conservatively keeping the error bars at 

the present level. 

Following (| Amati et al.ll2008h . a distribution mimicking 
the observed GRB redshift distribution in the range < 
z < 6 is used. A chosen percentage of these events was then 
replaced with redshifts uniformly distributed in the range 
6 < z < 10. This allows for a tuning of the number of 
high redshift GRBs in the mock sample that is used in the 
fits. Notice that different choices on the shape of the high 



redshift distribution of GRBs is approximately equivalent 
to a change on the redshift cutoff value and its percentage. 
Ultimately, only a significant number of high redshift GRBs 
can yield a sizable restriction on £o- The details about the 
actual low redshift distribution used to generate the mock 
data are discussed in the Appendix. 

Once a redshift distribution is obtained, lognormal dis- 
tributed values of E P) ; are attributed to each data point and 
are associated with a power law related Ei so - In this pa- 
per a value of 5.86 for log(E Pi i/lkeV) is used for the mean 
and a value of 1 is used for the variance. An extrinsic vari- 
ance (using (Text = 0.41) and Gaussian errors for E P; i and 
Eiso (20% for both) are then included, mimicking the current 
observational situation. It is verified that the results are not 
particularly sensitive to these choices. Several mocks with a 
varying number of low and high redshift GRBs where gener- 
ated and studied. The presented results consist of a typical 
mock sample with 500 events generated and 10% of high 
redshift GRBs, in a ACDM universe. 



4 COSMOLOGICAL DATA 

In order to gauge to which extent GRB data can constrain 
cosmological parameters, one confronts it with other well 
known cosmologically relevant observational tests such as 
SNe, BAO and the CMBR shift parameter. 



4.1 Supernovae 

The SNe sample from th e Sloan Digital Sk y Survey II 
(SDSS-II) first-year results ()Kessler et a.l.ll2009l ) is used, con- 
sisting of 288 SNe la with redshifts up to 1.55. SNe are 
used as distance indicators by comparing the theoretical dis- 
tance modulus, fith, obtained from the measured redshift in 
a given model of cosmological evolution, and with the in- 
ferred distance modulus, ^i bs, computed from fits to the 
SN e light curves (usin g the data of the mlsc2k2 fits found 
in (|Kessler et alj|2009h ). 

Specifically, the theoretical distance modulus is given 

by 



Atth = 51og£> L + 51og( — ) + 25 
Ha 



(10) 



with D L being the scaled (H independent) luminosity dis 
tance in Mpc, 

dz' 



D L {Z) = (1+2) 



E(z') 



(11) 



For each cosmological parameter choice, the used like- 
lihood is given by 



2 

XSN 



SN 
objects 



(^obs — fJ-th) 2 



(12) 



+ o"dis P + CT z i s the measurement variance 



where a 

for /.tobs, including the error from the fit, an intrinsic disper- 
sion error of <7dis P = 0.16, and a redshift error to account for 
the host galaxy peculiar movement and spectroscopic mea- 
surement. Only the SNe results are marginalized over Ho 
with a flat prior; all the other constraints assume a constant 
h = 0.7. 
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BAOs 



CMBs 



4.2 Baryon Acoustic Oscillations 

One also considers the constraints from the effect of the 
baryon acoustic peak of the large scale correlation function 
at lOO/i" 1 Mpc sepa ration detected by the SPSS Luminous 
Red Galaxy sample (|Reid et a l, 2009; Moldenhauer fc Ishakl 
2009). The peak position is related to the quantity 



A= yfih^E^)- 1 '' 6 



dz 



-i 2/3 



(13) 



zi Jo E(z)_ 

measured to be Ao — 0.493, with an error of a a = 0.017. 
The used likelihood is given by 



2 

Xbao 



Aq - A 



(14) 



4.3 



Cosmic Microwave Background Radiation 
shift parameter 



Here, the constraints f rom the CMBR WMAP7 observations 
Komatsu et al] 1201(f) are considered. The shift parameter 



Bond et al]|l997t r can be used as a distance prior to con- 



strain a given dark energy model. The shift parameter is 
given by the equation 



Rth 



(15) 



with Dl being the luminosity distance defined in Eq. 
and 2* the redshift at decoupling. The sta ndard fitting for- 
mula for z* is used (|Hu fc Sugivamalll996l ). This theoretical 
prediction is then constrained through the fitted WMAP7 
observation, i? bs = 1.725 with error ctr — 0.018, through 



2 _ / Rth — Ro 

XCMB — 

V OR 



(16) 



order to show what can be accomplished from a large pop- 
ulation of GRBs, despite the current level of measurement 
and theoretical uncertainties. A ACDM universe has been 
used to generate a mock sample (see section l3Tj) . This time, 
the results are marginalized over Q,u n with a flat prior in the 
interval 0.2 < fi Mo < 0.4. In Fig.lAlflthe results of the SNe 
and GRB constraints are shown. The constraints obtained 
are Co e [-0.21,0.82] and w £ [-0.84,-0.58], at 68% con- 
fidence level. From the SNe data alone, only a lower bound 
can be derived for Co in the considered parameter range. 

Once again, this SNe degeneracy in Co can be lifted 
by combining the SNe and CMBR constraints, yielding the 
bounds Co £ [-0.29, 0,181, ar^^ R £^-l. 08, -0.65] at 68% 
confidence level. Fig. IA3I combines all the available con- 
straints, namely GRB, SNe, BAO and CMBR. The CMBR 
data clearly provides the tighter constraints on all the pa- 
rameters. It is, however, highly degenerate in w and, hence, 
the SNe constraints on w are required to yield significant 
bounds. Notwithstanding, the GRB data has a similar pro- 
file to the CMBR bounds, and they can provide a signif- 
icant bound in Co- Thus, GRB data, in combination with 
SNe, provides an independent and compatible constraint on 
Co and w. On the other hand, the BAO result clearly does 
not yield a strong constraint on the value of either w or 
Co (as opposed to the SNe data), and does not significantly 
improve the results obtained from GRBs or CMBR. The 
overall combined results give the bounds Co £ [—0.10,0.08] 
and w £ [—0.89, —0.70] also at 68% confidence level. 

The GRB constraint on Co comes mostly from its high 
redshift valued data, whereas the low redshift valued SNe 
data presents a degeneracy in Co- Notice however that for 
the GRB and CMBR data, one encounters a degeneracy 
in Co in the form of a bend occurring at 77 < (i.e. Co > 
— 3w Q.de ). This occurs as for negative r\ and for high z, the 
evolution is dominated by the dark energy density, rendering 
the luminosity distance virtually independent of r\. 



5 RESULTS AND CONSTRAINTS 



One starts with the co nstraints obtained from the real ob- 
served 95 GRBs data (lAmati et aJ.ll200Sl. 120091) . combined 
with the SNe SDSS-II data (jKessler et all 120091 '). For this, 
O>m = Qdm + f2& is fixed at Qm = 0.3 with fij, = 
0.0445 and the Hubble parameter held at h = 0.7. The 
result can be seen in Fig. IA1I At 68% confidence level, 
the SNe data alone provides the limits Co £ [—1.93, 1.02] 
and w £ [—1.07, —0.62]. These results are, by themselves, 
highly degenerate in Co. Including the GRB real data, one 
improves slightly the constraints to Co £ [—1.15, 1.66] and 
w £ [-0.94,-0.58]. As expected, the status of the GRB 
data at present does not allow for a significant improve- 
ment over the SNe constraints on these two parameters 
(see Fig. [Al]). This SNe de generacy is usually lifted com- 
bining the SNe data with the CMBR WMAP7 shift param- 
eter observations yielding the bounds Co £ [—0.01, 0.13] and 
w £ [—0.83, —0.65] at 68% confidence level (not shown). 
Even though the present GRB data cannot compete with 
these CMBR constraints, with additional data they can pro- 
vide an important independent method of lifting the SNe 
degeneracy in Co- It must be stressed that these results are 
obtained for Qm fixed at 0.3. 

To illustrate this point, a mock population is chosen in 



6 DISCUSSION AND CONCLUSIONS 

In this work the use of GRB s as cosmological tools is 
considered. It has been shown (|Bertolami fc Tavares!l2006l : 
lAmati et al that GRBs have a great potential to mea- 

sure the value of Qdai independently from the CMBR con- 
straints. In the current work, it is shown that the present 
GRB data already gives a better constraint on the dark en- 
ergy and dark matter coupling parameter Co than the BAO 
results. The SNe results provide good bounds on w, but are 
more degenerate in Co- Despite of that, the experimental 
GRB sample is still too small to provide significant con- 
straints on Co- The combined result for SNe and real GRB 
data further illustrates this, yielding Co £ [—1.15, 1.66], with 
a width of ACo = 2.81, as opposed to the combined SNe 
and CMBR data limit Co £ [-0.01, 0.13] (ACo = 0.14), both 
at 68% confidence level. These results are obtained without 
marginalization, for fixed CIdm = 0.3. 

A mock population of GRBs was then generated show- 
ing that as the number of available events increases, GRB 
data becomes a more and more valuable tool in constraining 
the parameter Co - The GRBs complement the SNe constraint 
in a similar way the CMBR does. For a fixed fijvf = 0.3, 
the SNe and mock GRB data yields ACo = 0.82. 
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For a deeper insight on the future cosmological implica- 
tions of the GRB observations, the analysis was extended to 
a marginalization over £Im £ [0.2,0.4]. The combined SNe 
and CMBR results are ACo = 0.47 (68% CL), while for SNe 
and the mock GRB yields ACo — 1.03. Granting that the ob- 
tained bounds are not as accurate as the CMBR ones, they 
still provide a valuable independent measurement of these 
cosmological parameters. 

With the same marginalization in Qm , the combined 
result for SNe, CMBR and BAO is Co e [-0.27,0.13] (95% 
CL). The updated observations imp rove the previous re- 
sult Co G [-0.4,0.1] (|Guo et all 120071 ). Note, however, that 
tighter priors are used in this work and that this model is 
slightly different, with the inclusion of non-interacting bary- 
onic matter. 

It is interesting to compare the present results with 
the ones arising from estimates of the departure from the 
virial equilibrium of the Abell cluster A586. The bounds 
for r) from the Abell Cluster A586 yield rj G [ 3;.§§^dlQafeoc 
(Arj = 0-35), with w = -1 and z = 0.1708 and ft Mo = 0.28, 
|Bertolami et all I2007T I. Using SNe and mock GRB data 
from the present work and fixing Qm = 0.28, one ebssaiS-SOC 
ters tj 6 [0.93, 2.48] (A-q = 1.55) at 68% confidence level. 



If one considers the particular case of the GCG, the Abell 
Cluster A586 limits the a parameter to a G [0.21, 0.33] (68% 
CL). This compares with the combined SNe and mock GRB 
results, a e [0.25, 0.83] (68% CL). One should bear in mind 



that the GRB results were obtained from a mock population 
and are, thus, only indicative of what can be expected when 



the number of observed GRBs increases. The limit| dfeffpftme 



with the current observational SNe and CMBR data yield 
a e [0.03, 0.06 



I ti200 o I bento 2003 
In this work, the Amati correlation l)Amati et al.ll2008l ) 



has been used, given that it requires only two parame- 
ters whose determination can be inferred and increasing the 
number of useful GRB events available. Progress in the cali- 
bration and on the theoretical framework o:" this calibration 
would be invaluable, reducing the error margins in the GRB 
data and considerably improving the constraints on the pa- 
rameters. 
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APPENDIX A: GRB MOCK DATA 

For the mock sample generation the f ollowing d 
is considered (|Porciani fc MadadlioOlh 

dV(z) Rgrb(z] 





Hul996 


kamenshcik 



dN(Px < P s£ P 2 ) = 



dz 



l + z 



'Pi,* 



| SneSDSS 




wmap7 


1 P 


^zhang 




Drimer 



(Al 



where dV/dz is the comoving volume element, Rgrb is the 
comoving GRB rate density and e(P) is the detector effi- 
ciency as a function of photon flux. The quantity ifi(L) ipEMDpri 



normalized GRB luminosity function and L is a "isotropic 
equivalent" burst luminosity L = j^°^ e y eV E S(E) dE, for 
the energy E and S(E) is the rest-frame photon luminosity 
of the source. 

The generation process is performed as follows: using a 
Monte-Carlo generator, a sample of the desired number of 
GRBs with z in the range < z < 6 is generated, using 
Eq. (|A1[) . Then, a desired percentage is replaced randomly 
by events in the range 6 < z < 10, using a flat distribution. 
The Uiso is randomly attributed according to a Gaussian dis- 
tribution and _E Pi i is then calculated using a power law with 
the parameters calculated using a fit of the real GRB sam- 
ple. The errors are added afterwards assuming a Gaussian 
distribution. 
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Figure Al. Constraints on the values of w and £o obtained from 
the SNc data (dashed line) and its combination with the 95 GRB 
data (full line). The 68% and 90% confidence levels are depicted. 
One fixes: SImo = 0.3 and h = 0.7. 




Figure A3. The combined constraints on the values of w and £o 
obtained from the SNe (red), BAO (light blue), CMB (blue) and 
a mock 500 GRB population (green). The GRB and SNe contours 
are the same as in Fig, IA2l The CMB contour is the narrower one 
(in blue) and the BAO contour is in background (light blue) . The 
dark patch shows the combined region. The 68% and 95% confi- 
dence levels are shown. A marginalization over Qm € [0.2,0.4] 
has been carried out. 



fig:SNBA0GRBCM 



Figure A2. Constraints on the values of w and £o obtained from 
SNe and a mock of 500 GRB population. The wider contour in 
w corresponds to the GRBs (in green, for online readers) and 
the wider contour in to the SNe (in red, for online readers). 
Both 68% and 95% confidence levels are shown. The shaded area 
refers to the combined confidence region. A marginalization over 
SIm G [0.2,0.4] has been considered. 
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